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ABSTRACT 

Comparing the gravitational acceleration induced on the Local Group of 
galaxies by different tracers of the underline density field we estimate, within 
the linear gravitational instability theory and the linear biasing ansatz, their 
relative bias factors. Using optical SSRS2 galaxies, IRAS (PSCz) galaxies and 
Abell/ACO clusters, we find b ,i « 1.21 ± 0.06 and b c ,i ~ 4.3 ± 0.8, in 
agreement with other recent studies. Finally, there is an excellent one-to-one 
correspondence of the PSCz and Abell/ACO cluster dipole profiles, once the 
latter is rescaled by bc,i, out to at least ~ 150 h~ x Mpc. 
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1 INTRODUCTION 

Different classes of extragalactic objects trace the un- 
derlying matter distribution differently. The realization 
of such a behaviour arose from the fact that the am- 
plitude of the 2-point correlation function of clusters of 
galaxies is significantly higher than that of galaxies (cf. 
Bahcall & Soneira 1983). This was suggested by Kaiser 
(1984) as a result of the clustering characteristics of dif- 
ferent height peaks in an underlying random Gaussian 
field. A first order description of the effect is provided 
by linear biasing in which the extragalactic mass tracer 
fluctuation field is related to that of the underlying mass 
by: 

(<WP)traccr = b (<Wp)mass W 

with b the linear bias factor. Even in this simplistic 
model the bias cannot be measured directly and only 
theoretical considerations or numerical simulations can 
provide some clues regarding its value. However, the rel- 
ative bias between different tracers can be measured 
and such attempts, using their clustering properties, 



clusters: general - large-scale structure of universe - in- 

have provided interesting, although somewhat conflict- 
ing, results. Lahav, Nemiroff & Piran (1990) compar- 
ing the angular correlation function of different sub- 
samples of the UGC, ESO and IRAS catalogues find 
an optical to IR galaxy bias factor, bo,i, ranging from 
1 to 2 with preferred value bo,i ~ 1.7. Babul & Post- 
man (1990) using the spatial correlation function of the 
CfA and IRAS galaxies find bo,i — 1.2, while compar- 
ing the QDOT correlation function (Saunders, Rowan- 
Robinson & Lawrence 1992) with that of APM galaxies 
(Maddox et al. 1990) one finds bo,i ~ 1.4. Similarly, 
Oliver et al. (1996) comparing the clustering properties 
of the APM-Stromolo survey of optical galaxies and an 
extended IRAS redshift survey found bo.i ~ L2 ± 0.05. 
Strauss et al. (1992a) using the 1.936 Jy IRAS sample 
find that the overdensity ratio between CfA and IRAS 
galaxies within a sphere centered on Virgo with the Lo- 
cal Group on the periphery gives bo ,i — 1.2 while their 
correlation function analysis provides discrepant results 
when comparing IRAS to CfA or SSRS optical galaxies 
(with bo, i ~ 2 and 1, respectively). Recently, Willmer, 
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daCosta & Pellegrini (1998) using the SSRS2 sample 
of optical galaxies and comparing with the clustering 
properties of the 1.2 Jy IRAS survey find bo,i ~ 1.2 
and 1.4 (±0.07) in redshift and real space respectively, 
while Seaborne et al (1999) performing a similar analysis 
between the PSC and Stromolo-APM redshift surveys 
find bo, i ~ 1.3 ± 0.1. 

A different approach using the dynamics of the lo- 
cal group of galaxies, was proposed in Plionis (1995) 
and Kolokotronis et al. (1996). Traditionally, dynamical 
studies have been used in an attempt to constrain the 
value of the cosmic density parameter, Q a , by assuming 
linear theory and comparing observed galaxy or clus- 
ter peculiar velocities with estimated accelerations (cf. 
Strauss & Willick 1995). However due to biasing only 
the combination Q,®' 6 jb can be estimated. Such an anal- 
ysis has been extensively applied to the Local Group 
of galaxies, since its peculiar velocity is accurately de- 
termined from the CMB temperature dipole (Kogut et 
al. 1996) and its gravitational acceleration can be mea- 
sured from the dipole moment of the surrounding spa- 
tial distribution of different mass tracers. Within lin- 
ear theory acceleration and peculiar velocity should be 
aligned and this indeed has been found to be the case 
using optical, IR galaxies, X-ray or optical cluster sur- 
veys and AGN's (cf. reviews of Strauss & Willick 1995, 
Dekel 1997 and references therein). In the linear biasing 
framework the different mass tracers should therefore 
exibit similar dipole profiles differing only in their am- 
plitudes, the ratio of which is a measure of their relative 
bias. Therefore, one can estimate the relative bias fac- 
tor between different mass tracers, because in the inter- 
comparison of their velocity-acceleration relations the 
Qo parameter as well as the velocity cancels out. 

In this study we use the recently completed PSCz 
IRAS galaxy survey (Saunders et al. 1999), the SSRS2 
optical galaxy catalogue (DaCosta et al. 1998) and a 
subsample of the Abell/ACO cluster catalogue (as de- 
fined in Branchini & Plionis 1996) to estimate their rel- 
ative bias factors by comparing their dipole moments. 

2 METHOD 

Using linear perturbation theory one can relate the 
gravitational acceleration of an observer, induced by the 
surrounding mass distribution, to her/his peculiar ve- 
locity: 

n - 6 If r tt os 

^ = ^lj s ^W dr = -ir D{r) (2) 

The dipole moment, D, is estimated by weighting the 
unit directional vector pointing to the position of each 



tracer, with its gravitational weight and summing over 
the tracer distribution; 

D 47r(n) (j>(r) r 2 * ^ 
with 

4>{r) = / '" aX *(L) dL (4) 

where $(1/) is the luminosity function of the objects un- 
der study, Lmi n (r) = 47rr 2 5'ii m , with Sum the flux limit 
of the sample and (n) is the mean tracer number den- 
sity, given by integrating the luminosity function over 
the whole luminosity range. 

Using two different tracers, i and j, of the underly- 
ing matter density field to determine the Local Group 
acceleration one can write: v(r) = Q°- 6 Di(r)/bi = 
Q.°- 6 r>j(r)/bj and therefore we can obtain an estimate 
of their relative bias factor from: 

Mr) = |W = §tW ( 5 ) 

Since the dipole is a cumulative quantity and at each 
distance it depends on all previous shells, we cannot 
define an unbiassed x 2 statistic to fit eq|]. Rather, we 
can obtain a crude estimate of the reliability of the re- 
sulting bias factor by estimating Pearson's correlation 
coefficient, Rij, between the two dipole profiles (see 
Kolokotronis et al. 1996); a value Rij ~ 1 would in- 
dicate a perfect match of the two dipole profiles and 
thus a very reliable estimate of their relative linear bias 
factor. 

A statistically more reliable approach is to assume 
that the differential dipoles, estimated in equal volume 
shells, are independent of each other and then fit bij 
according to: 

2 _ \ ^ [Dj,h — bijDj.k — Ck) , . 

X ~ 2—t CT 2 + b 2 a 2 ^ ' 

where C is the zero-point offset of the relation and a 
is the corresponding shot-noise errors, estimated by us- 
ing either of two methods; a Monte-Carlo approach in 
which the angular coordinates of all tracers are random- 
ized while keeping their distance, and thus their selec- 
tion function, unchanged or the analytic estimation of 
Strauss et al. (1992b); a 2 ~ J2 ^"V" 4 ^- 1 + 1). 

3 DATA 

We use in our analysis three different catalogues of mass 
tracers; 

• The recently completed IRAS flux-limited 60-^m 
redshift survey (PSCz) which is described in Saunders 
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et al. (1999). It is based on the IRAS Point Source Cat- 
alogue and contains ~ 15000 galaxies with flux > 0.6 
Jy. The subsample we use, defined by |6| > 10° and lim- 
iting galaxy distance of 180 ft -1 Mpc, contains ~ 10097 
galaxies and covers ~ 82% of the sky. 

• The SSRS2 catalogue of optical galaxies (DaCosta 
et al. 1998) which is magnitude limited to tub = 15.5 
and contains 3573 galaxies in the South (—40° < S < 
-2.5°, b < -40°) and 1939 galaxies in the North (S < 
0°, b > 35°), covering in total 13.5% of the sky. 

• A volume limited subsample of the Abell/ACO 
cluster catalogue, with |6| > 10° and limited within 180 
h^ 1 Mpc (see Branchini & Plionis 1996). Our sample 
contains 197 clusters. 

3.1 Determining distances from redshifts 

All heliocentric redshifts are first transformed to the 
Local Group frame using cz ~ czq + 300sin(Z) cos(6). 
We then derive the distance of each tracer by using: 

r = -jj£-(l-(l + «- Sz)- 1/2 ) (l + z- Szf /2 (7) 
Ho 

where H = 100 h Mpc and 8z is a non-linear term to 
correct the redshifts for the tracer peculiar velocities: 

5z=-(u(r)-u(0))-r (8) 

c 

with u(0) the peculiar velocity of the Local Group and 
u(r) the peculiar velocity of a galaxy or cluster at po- 
sition r. Instead of using elaborate 3D reconstruction 
schemes (cf. Schmoldt et al 1999; Branchini & Plionis 
1996; Branchini et al 1999; Rowan- Robinson et al. 1999) 
to estimate this term, we decided to use a rather simplis- 
tic velocity field model (see Basilakos & Plionis 1998) 
to treat consistently all three data sets (a self-consistent 
3D reconstruction of the SSRS2 density field is in any 
case not possible due to the small area covered by the 
survey). Our simplistic velocity field model was found 
in Basilakos & Plionis (1998) to be sufficient in order 
to recover the IRAS 1.2Jy and QDOT 3-D dipole. We 
remind the reader the main assumptions of this model: 

(a) The tracer peculiar velocities can be split in two 
vector components; that of a bulk flow and of a local 
non-linear term: 

u(r) = V 6 „ !fc (r)+u n i(r) (9) 

(b) The first component dominates and thus that 

u(r)-f «V 6 nifc(r)-f (10) 

We then use the observed bulk flow direction and pro- 
file, as a function of distance, given by Dekel (1997) 
and combined with that of Branchini, Plionis & Sciama 
(1996). The zero-point, Vbuik{0), and the direction of the 
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Figure 1. Mean SSRS2 (open symbols) and PSCz (filled 
symbols) galaxy space density estimated in equal volume 
shells with its Poissonian uncertainty. The continuous line 
is the SSRS2 density reduced by a factor 2. Based on 13.5% 
of the sky (SSRS2 survey area). 

bulk flow is estimated applying eq.([]) at r = and as- 
suming, due to the "coldness" of the local velocity field 
(cf. Peebles 1988), that u n i(0) ~ u inf = 200 km/sec 
(where Ui n t is the LG infall velocity to the Virgo Super- 
cluster) . 

3.2 Galaxy densities 

To estimate the local acceleration field it is necessary 
to recover the true galaxy density field from the ob- 
served flux-limited samples. This is done by weighting 
each galaxy by 4>~ 1 (r) , where (f>(r), is defined in eq^. For 
the PSCz sample we use the Saunders et al. (1990) lumi- 
nosity function derived from the QDOT catalogue, with 
= 7.5 x 10 7 h~ 2 L since lower luminosity galax- 
ies are not represented well in the available samples (cf. 
Rowan- Robinson et al. 1990), and L m ax = 10 13 K~ 2 Lq. 
For the SSRS2 sample we use the Schechter luminosity 
function of Marzke et al. (1998) with M max — —22 and 
M min — —13.8. 

In figure 1 we present the mean density and its 
Poissonian uncertainty of PSCz and SSRS2 galaxies in 
their common area (that of the SSRS2 sample) and in 
equal volume shells (with SV ~ 4.5 x 10 6 ' h~ 3 Mpc 3 ). 
Their densities are extremely comparable, differing only 
by a constant factor ({po) / (pi) = 2.03 ± 0.16). 

4 RESULTS 

4.1 Optical to IR galaxy bias 

We first present the results of the intercomparison of the 
SSRS2 and PSCz samples in their common sky area and 
for r > 15 h" 1 Mpc. Infi gure 2a we show the amplitudes 
of the two dipoles as a function of distance from the 
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LG. The monotonic dipole increase reflects the fact that 
we are measuring only the component of the whole sky 
dipole which is due to the particular area covered by 
the sky restricted SSRS2 sample. It is apparent that 
the shapes of the two dipole amplitudes are extremely 
similar, giving correlation coefficient R ~ 0.97. 

In figure 2b we present the direct dipole ratio (eq.^) 
in the LG frame (open symbols), while as filled squares 
the results of the fit of eq.[], as a function of maximum 
distance used. No significant differences are found when 
correcting distances for peculiar velocities. It is evident 
that the different estimates are consistent with each 
other, especially for r > 50 h~ x Mpc where the direct 
dipole ratio becomes flat. It is essential, however, to ver- 
ify whether such a good dipole-profile correlation could 
result solely due to the small solid angle used, ie., to 
investigate whether the survey geometry, coupled with 
the galaxy selection function, dominates the dipole sig- 
nal of both samples. To this end we have generated 100 
mock SSRS2 samples by reshuffling the galaxy angu- 
lar coordinates while leaving their distances, and thus 
selection function, unchanged. If the reshuffled dipole 
profile resembles that of the SSRS2 original one, then 
this would indicate the existence of the previously sug- 
gested bias. In figure 3 we present both (a) the dipole 
ratio between the original SSRS2 dipole and the reshuf- 
fled one together with its scatter and (b) the SSRS2 and 
PSCz dipole ratio with the latter rescaled by boj. It is 
evident that the former ratio deviates significantly from 
one, an indication that our comparison is not dominated 
by the suspected bias. 

Using the differential equal volume dipoles to fit 
eq.| for 10 < r < 185 h' 1 Mpc, we find b ,i ^ 
1.24 ± 0.04, with zero-point C ~ -105 ± 50 km/sec and 
\ 2 ~ 6.3 for 6 degrees of freedom. The fit is performed 
using the FITEXY routine of Press et al. (1992) and the 
uncertainties of the fitted parameters correspond to the 
A\ 2 = 1 confidence region boundary. The small zero- 
point offset could be due to uncertainties in tracing the 
very local contributions to the LG dipole. Indeed, in- 
tegrating the the SSRS2 and PSCz dipoles for r > 15 
h~ x Mpc (presented in figure 2) we find no zero-point 
offset C ~ —40 ± 45 km/sec but a slightly smaller bias 
factor b j ^ 1.17 ± 0.04, with \ 2 ~ 7.5 for 6 degrees 
of freedom. We derive a mean estimate of the bias fac- 
tor and its uncertainty by varying both the inner and 
outer dipole integration limits in eq.^ and by taking 
into account the slight differences of the results based 
on the differential dipole (eq.|E]) . The resulting optical to 
IR galaxy bias factor is: 
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Figure 2. (a) Dipole amplitude comparison for the SSRS2 
and PSCz samples, (b) The estimated bias from comparing 
cumulative dipoles using eq|j (open circles) and differential 
dipoles using eq.|(j (filled squares). 
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Figure 3. Cumulative Dipole Ratio test: Solid line repre- 
sents the ratio of the original SSRS2 and mean reshuffled 
dipoles; with its scatter estimated over 100 reshufflings of 
the galaxies angular coordinates. Solid points represent the 
ratio of the SSRS2 and PSCz dipoles with the latter rescaled 
by bo,i — 1.17 (see text) while the errors represents the 
propagation of the shot-noise dipoles. 



Our result is in quite good agreement (within la) with 
that of Seaborne et al (1999), which is based on the 
PSCz and APM galaxy clustering properties on rela- 
tively smaller scales than those probed by our dipole 
analysis. 

4.2 Rich cluster to IR galaxy bias 

In the case of the cluster and PSCz samples we have a 
nearly full sky coverage, except at low-galactic latitudes. 
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In order to recover the whole sky dipole we use a spher- 
ical harmonic approach to "fill" the unobserved part of 
the sky (cf. Lahav 1987; Plionis & Valdarnini 1991). 
This approach has been found to provide compatible 
results with the cloning and interpolating method (see 
Branchini & Plionis 1996 for such a comparison in the 
context of the cluster dipole). 

The main drawback of comparing the cluster and 
galaxy dipoles arises from the fact that the Abell/ACO 
cluster distribution is incomplete in the local universe, 
since it does not include the Virgo cluster, an important 
contributor of the local velocity field (cf. Tully & Shaya 
1984). Therefore the direct comparison of the dipole am- 
plitudes is hampered by this zero-point uncertainty. We 
can attempt to correct for this problem by: 

• including the local Virgo contribution to the cluster 
dipole by assigning an appropriate Abell number count 
(Na) weight to the Virgo cluster, 

• excluding from the PSCz dipole the very near contri- 
butions (< 8 h' 1 Mpc). 

A first attempt to derive the cluster to IR galaxy bias, 
using such a procedure and the Abell/ACO and QDOT 
catalogues, was present in Plionis (1995). If clusters and 
galaxies do trace the same underline field, as indeed ap- 
pears to be the case (cf. Branchini et al. 1999; Bran- 
chini, Zehavi, Plionis & Dekel 1999), then we should be 
able to fit the two profiles, using eq|(| varying the Virgo 
cluster weight. The appropriate value of Na is that for 
which the zero-point, C, of the fit vanishes and thus we 
will consider as our preferred bias parameter the cor- 
responding value of be, i- Statistically, this procedure 
does not provide a rigorous significance indication, due 
to the fact that the dipoles are cumulative quantities, 
but it provides a means of comparing quantitatively the 
two dipole profiles. In order to test the robustness of the 
resulting bias parameter we fit the two dipole profiles as 
a function of distance. 

In figure 4a we present the resulting bias parameter 
versus the zero point, C, for the Virgo cluster weights 
that provide a fit with C ~ 0. The different connected 
point arrays correspond to results based on the different 
Na weights while different points in each array corre- 
spond to different upper distance limits used for the fit, 
which increase as indicated by the arrow. Taking into 
account that the zero-point uncertainty is SC ~ 130 
km/sec we conclude that the Virgo cluster weight for 
which C ~ is Na ~ 24 ± 4, confirming the notion 
that Virgo corresponds to a richness class R — Abell 
cluster. A consistency check, that we pass with success, 
is that the Virgocentric infall velocity that corresponds 
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Figure 4. (a) Fitted bias parameter between clusters and 
IRAS galaxies versus zero-point, C, for different Virgo cluster 
weights. The different points within each connected array 
represent the results of the fit to different limiting distance 
(increasing in the direction of the arrow), (b) The fitted bias 
parameter (cq.^j) as a function of limiting distance for N a = 
24 (continuous line). The broken line is the corresponding 
direct dipole ratio (eq.tj). 



to this weight is Mi n f — 300 ± 40 km/sec, a value in 
agreement with most available determinations. 

In figure 4b we present the fitted bias parameter as 
a function of upper distance limit (points) and the direct 
dipole ratio (eq.^) as broken lines for the Na = 24 case. 
Both seem consistent within their uncertainties, espe- 
cially for distances ~ 140 — 150 h" 1 Mpc, which roughly 
corresponds to the apparent cluster dipole convergence 
depth. The main result regarding the bias parameter is: 

b c ,i ~ 4.3 ±0.8 , 

which interestingly is mostly independently of the Virgo 
cluster weights (as can be seen in figure 4a), since such 
differences are absorbed in the value of C. The uncer- 
tainty in be, i reflects (a) variations due to different 
Virgo weights, (b) the variation between the eq.^ and 
eq.[i| solutions and (c) the scatter around these solutions 
(see figure 4b). It does not include, however, the cosmic 
variance which results from the use of only one observer. 
Our results are in excellent agreement with those of Pea- 
cock & Dodds (1994) and Branchini et al. (1999b) based 
on completely different approaches. 

In figure 5 we present a direct comparison of the 
PSCz and Abell/ACO cluster dipoles, out to 150 h' 1 
Mpc, after having scaled down the latter by bc,i = 4.3. 
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Figure 5. Comparison between the PSCz and the 
ABELL/ACO dipoles, after scaling down the latter by a bias 
factor of 4.3. 



The errorbars in the cluster dipole reflect mainly the un- 
certainty of the cluster density between the Abell and 
ACO parts of the sample (see Plionis & Kolokotronis 
1998 and references therein). The two profiles are in ex- 
cellent agreement, at least, up to ~ 150 hT x Mpc with 
correlation coefficient R — 0.86. This is a further in- 
dication that the two density fields are consistent with 
each other out to these distances and supports the exis- 
tence of dipole contributions from large depths (see also 
Schmoldt et al. 1999; Basilakos & Plionis 1998), sug- 
gestions which were first put forward by Plionis (1988), 
on the basis of the Lick counts, by Rowan-Robinson 
et al. (1990) on the basis of the QDOT survey and by 
Scaramella et al. (1991) and Plionis & Valdarnini (1991) 
on the basis of Abell/ACO clusters. A thorough inves- 
tigation of the deep PSCz dipole (distances > 150 hT 1 
Mpc) will be presented in Rowan- Robinson et al. (2000). 



5 CONCLUSIONS 

We have used a novel approach, based on the Local 
Group dipole properties, to estimate the relative bias 
parameter of different mass tracers. We find that the 
optical to IR galaxy bias parameter is bo,i — 1.21±0.06, 
while the rich cluster to IR galaxy bias is bc,i — 4.3±0.8. 
Our results are in good agreement with others based 
on different approaches. We find that the IR galaxy 
and rich cluster dipole profiles are extremely compat- 
ible once the latter is rescaled by bc,i out to at least 
~ 150 ft" 1 Mpc. 
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